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ABSTRACT 

We consider the effects of self-gravity on the hydrostatic balance in the vertical direction of a gaseous 
disk and discuss the possible signature of the self-gravity that may be captured by the direct imaging 
observations of protoplanetary disks in future. In this paper, we consider a vertically isothermal disk 
in order to isolate the effects of self-gravity. The specific disk model we consider in this paper is the one 
with a radial surface density gap, at which the Toomre's Q-parameter of the disk varies rapidly in the 
radial direction. We calculate the vertical structure of the disk including the effects of self-gravity. We 
then calculate the scattered light and the dust thermal emission. We hnd that if the disk is massive 
enough and the effects of self-gravity come into play, a weak bump-like structure at the gap edge 
appears in the near-infrared (NIR) scattered light, while no such bump-like structure is seen in the 
sub-mm dust continuum image. The appearance of the bump is caused by the variation of the height 
of the surface in the NIR wavelength. If such bump-like feature is detected in future direct imaging 
observations, with the combination of sub-mm observations, it will bring us useful information about 
the physical states of the disk. 
Subject headings: protoplanetary disks 



1. INTRODUCTION 

Direct imaging observations of protoplanetary disks 
have been making rapid progress. In the near-infrared 
(NIR) wavelength, the strategic observation program 
using Subaru telescope, Subaru Strategic Exploration 
of Exoplanets and Disks (SEEDS), is now under way 
(Tamura 2009). Several results to date reveal that disks 
have diverse structure s such as an inner hole, a gap or a 
ring (jThalmann et all [20101: Hashimoto et all 120111 ). In 
the sub-millimeter (sub-mm) wavelength, Atacama Large 
Millimeter/submillimeter Array (ALMA) is expected to 
observe disks at a spatial resolution that is comparable 
to, or better than, NIR wavelength. 

Direct imaging observation is useful in probing the de- 
tailed structures in the outer disk. For example, a tele- 
scope with 0.1 arc-second resolution has an ability to 
resolve the structure with ~ 10AU scale if the disk is 
located at lOOpc. Such a spatial resolution is interesting 
since 10AU is comparable with the disk scale height at 
100AU away from the central star. The disk scale height 
is an important scale in considering the dynamical pro- 
cesses of disks. 

From a theoretical point of view, it is important to 
make a disk model that can be used for the interpre- 
tation of direct imaging observations of protoplanetary 
disks, and discuss what we can derive from the observed 
structures. Comparing the information from different 
wavelengths is essentially important in order to reveal 
the real structure of the observed disks. 

In this paper, we investigate the effects of self-gravity 
on the disk vertical structure in the hydrostatic equilib- 
rium and discuss how the signature of self-gravity ap- 
pears in the direct imaging observations. We consider 
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the observations of a disk at NIR wavelength, which 
mainly observe the scattered light from the upper layer of 
the disk, and the sub-mm dust continuum observations, 
which reveal the surface density structure of the disk. 

One of the fundamental and important properties of 
the protoplanetary disk is its mass. The sub-mm dust 
continuum observation is widely used to estimate the 
mass of the disk, although it suffers the ambiguity of 
the dust opacity, which results in at least a factor of two 
uncertainty in the disk mass. In this paper, we consider 
the effects of self-gravity on the hydrostatic balance of 
the disk vertical structure, and discuss how the effects 
of self-gravity may appear in the direct imaging observa- 
tions. If there is a feature that reflects the self-gravity, 
it would give an alternative implication about the disk 
mass other than the strength of the dust continuum. 

In order to isolate the effects of self-gravity in the disk 
structure, we consider a disk with a simple radial surface 
density variation, namely, a power-law profile with an 
axisymmetric gap. The variation of the surface density 
may be produced by the gravitational perturbation by 
an embedded planet, turbulence, or other causes, but we 
do not query the cause of such a structure. The question 
we address is that how such surface density variations 
may be observed by direct imaging observations at NIR 
or sub-mm wavelengths and we look for a feature that is 
unique to a self-gravitating disk. 

There are several works on the modeling of a self- 
gravitating disk for direct imaging observations, espe- 
cially in relation to the gravitoturbulence and t he re- 
sulting planet formation. iNarayanan et all (|2006f ) have 
performed a detailed calculations of line emission based 
on the h ydro d y nami c simulations of gravitational insta- 
bility by IBossI ()2001l ) , and have shown that unique fea- 
tures in HC O + line may appear clo s e to a forming gi- 
ant planet. Uang-Condell and Boss! (|2007l ) have calcu- 
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lated the scattered light using the same re s ults o f hydro- 
dynamic calculations as iNaravanan et all (|2006t ). They 
have shown that the vertically integrated surface den- 
sity and the scattered light have no correlation, since 
strong turbulence induced by a gravitational instability 
produces blobs in the upper layer of the disk, and the 
scattered light is strongly affected by them. 

In this paper, we consider a laminar disk where the 
vertical hydrostatic balance is achieved. Therefore, our 
model is not useful in a strongly turbulent disk where 
the Toomre's Q-parameter (defined in the later section) 
becomes smaller than ~ 1.5 — 2, when non-axisymmetric 
gravitational in stability cau s es the disk to be in a turb u- 
lent state (e.g., iMeiia et al.l (|2005| ) : iBolev et all p006)). 

The paper is constructed as follows. In Section [2] we 
describe the disk model and the methods of calculations. 
In Section [3J we show the results of calculations. In 
this section, we also describe the features that appear 
in the NIR scattered light image and are unique to self- 
gravitating disks. In Section [4] we discuss the observa- 
tional implications. Section[5]is for summary and further 



discussions. 



2. DISK MODEL 



In this section, we describe the basic setup of our 
model. We consider a model where the outer disk has 
a gapped structure. We consider a cylindrical coordi- 
nate system (r, <fr, z) with the origin at the central star. 
As a simple model that features the effects of self-gravity, 
we assume an axisymmetric disk with a gap: 



hance/gapV ) ^Yqq^jJ g/ cm j (1) 

where /enhance is a parameter that controls the over- 
all surface density of the disk and f gap (r) is a function 
that determines the shape of the gap. We vary the disk 
mass by changing /enhance from 1 to 8.15, but in the 
subsequent section, we especially consider the case of 
/enhance = 1 and 6.6 as illustrative examples. The case 
of /enhance = 6.6 is rather extreme, since the typical Q- 
value in this case is close to 2, where turbulence may take 
place. However, the effects of self-gravity on the vertical 
structure is prominent in this extreme case. We assume 
a Gaussian gap profile, 



/gap(r) = 1 - aexp 



A gap 



(2) 



with a = 0.95, r gap = 100AU, and A gap = 30AU. The 
parameters a, r gap , and A gap parameterize the depth, 
location and width of the gap, respectively. 

In this paper, in order to make the problem simple and 
tractable, we make a simple assumption for the temper- 
ature of the disk. We assume a vertically isothermal disk 
where the temperature T(r) is given by 



f r \ ~V 2 
r W = 30 (l00Au) K - 



(3) 



Although we do not question ourselves about the origin 
of such surface density and temperature structure in de- 
tail, we comment here that it is reported that a massive 
planet emedded in a disk causes a deep gap in the sur- 
face density profile, while the temperature does not vary 



as si gnificantly as the surface density (|D'Angelo et al.l 
[2001 . 

By varying the overall surface density by /enhance while 
keeping the temperature profile fixed, we vary the disk 
Toomre's Q-parameter, Q = Qc/irGS. The Q-parameter 
is also varied in the radial direction due to the gap. Fig- 
ure[T]shows one example of the disk model. In this figure, 
we show the profile of surface density and Q-value with 

/enhance 1 ■ 

From the surface density and the temperature profile, 
we obtain the vertical structure of the disk. We assume 
that the hydrostatic equilibrium is reached and that the 
disk is geometrically thin. We solve 



dp(z) 
dz 



p{z)l&z+^ 



dz 



(4) 



where p — c 2 p is the pressure and Q(r) = (GM^/r 3 ) 1 / 2 is 
the Keplerian angular frequency. The sound speed c(r) is 
obtained from the temperature profile given by equation 
© by 

cw = io5 (Ssr™ /8 ' < 5 » 

In equation denotes the gravitational potential of 

the disk, which is determined by Poisson equation 



V 2 i/>d = 4?rGp. 



(6) 



We solve equations ^ and (JSJ iteratively to find a con- 
sistent solution. If the variation of the surface density 
is slow in the r adial direction, the local approach by 
Paczynski (|1978l ) yields a reasonable model. However, 
in the model presented in this paper, such variation may 
be rapid due to the presence of the gap, so it is necessary 
to solve full Poisson equation. 

Once we obtain the density structure of the disk, we 
solve the equations of radiative transfer to obtain the 
observed image. We consider the thermal emission of the 
disk and the light from the central star scattered by the 
disk surface. The former is important in the observations 
in sub-mm wavelength and the latter is important in the 
NIR imaging observations. For simplicity, we assume 
that the dust and gas are well-mixed and the optical 
properties of the dust particles do not vary in the disk. 
We also assume that the disk is face-on. 

We follow the approach by [DAlcssio et al. (199^) in 
calculating the thermal emission and the scattered light. 
The thermal emission may be calculated by integrating 



dl 



the 



dZ 



= -npB(T)exp[-T{Z)\ 



(7) 



where Z is the coordinate along the line of sight, k is the 
absorption cross section per unit gas mass, B(T) is the 
Planck function, and r is the optical depth along the line 
of sight, which is calculated by 



dr 

dz = -™ 



(8) 



where x is the total extinction (including absorption and 
scattering) cross section per unit gas mass. In calculat- 
ing the scattered light, we assume a single and isotropic 
scattering. Then, the scattered light can be calculated 
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by integrating 

Jjscat 

= -aW{r, z)B{T*) exp [-r rad - r{Z)\ , (9) 

along the line of sight. Here, a = (x — k)p is the scat- 
tering coefficient, W(r, z) = l/4(r 2 + z' 2 ) (r is the radial 
coordinate and z is the vertical coordinate) is the geomet- 
ric dilution factor fe.g.. lJang-Condellll2009l ). B(T*) is the 
blackbody radiation of the central star, and r ra d is the 
optical depth for the stellar irradiation. In determining 
W[r), we assume that the central star is a point source, 
which is a reasonable assumption when considering the 
outer disk. 

The opacity of dust particles is one of the major uncer- 
tainties in the observational modeling of the disk. There- 
fore, we use a range of values of the opacity in order to 
find an opacity-independent signature. We use the total 
extinction x with 40 — 200cm 2 /g for the modeling of the 
NIR observation, and 1.2 x 10~ 4 — 6 x 10~ 3 cm 2 /g for the 
modeling of sub-mm observation. In the subsequent sec- 
tion, where we describe the effects of self-gravity in detail, 
we show the results of x = 100cm 2 /g for NIR observa- 
tions and x = 0.003cm 2 /g for sub-mm observations. We 
use the fixed values of scattering albedo, u> = {x~ K )/x 
for simplicity. We use uj = 0.9 for the NIR modeling and 
oj = 10~ 3 for the sub-mm modeling. 

3. SIGNATURE OF A SELF-GRAVITATING DISK 

Since we assume that the disk structure is axisymmct- 
ric and the disk is face-on, the resulting image is also 
axisymmetric. Figure [2] shows the examples of the ob- 
served radial profile of the relative intensity in NIR wave- 
length and sub-mm wavelength. The disk brightness is 
normalized by that at 50AU. The disk with / cn hancc = 1 
(low-mass disk) and the disk with /enhance = 6.6 (high- 
mass disk) are shown. The total extinction coefficient x 
of the disk is 100cm 2 /g in NIR and 0.003cm 2 /g in sub- 
mm, respectively. 

In the sub-mm observation, there is no significant dif- 
ference in the profile of the disk brightness between the 
low-mass disk case and the high-mass disk case. This is 
because in the sub-mm wavelength, the disk is optically 
thin, and the observed flux is proportional to the disk 
surface density. In our model, the profile of the surface 
density is the same in the low-mass and the high-mass 
case except for the factor /enhance- Therefore, the relative 
intensity profiles are essentially the same in both cases. 

In the NIR observation, on the other hand, we see a 
difference in the radial structure of the relative intensity 
between the low-mass disk case and the high-mass disk 
case. In the case of the low-mass disk, the profile of 
the relative intensity is very similar to that of sub-mm. 
There is a dark region in the disk around the place where 
the gap is present. In the case of the high-mass disk, we 
see that there is a bump in the relative intensity profile 
at ~ 80AU region (indicated by an arrow) beside the 
dark region around ~ 100AU. We also note that there is 
another bump-like structure outside the gap region. In 
Figure [21 we see such a structure around 120 — 140AU 
region. In the region further away from the central star, 
we see that the NIR brightness is weakened. 

We now focus on the bump-like structure that appears 
in the NIR profile just inside the gap region (indicated 



by an arrow in Figure [5]) , and we shall show that this 
structure indeed arises from the effects of self-gravity. 

In the left panel of Figure [3l we compare the results 
of NIR scattered light with and without self-gravity. We 
use the high-mass disk model with /enhance = 6.6. The 
solid line shows the same calculation with the dashed 
line in the left panel of Figure [5] and the dashed line in 
Figure [3] shows the results for the same disk model but 
the self-gravity is artificially turned off. We see the bump 
structure around ~ 80AU region in the calculations with 
self-gravity but it is not significant in the calculations 
without self-gravity. 

In order for the bump structure to be seen, it is also 
necessary to have a gap in the surface density structure. 
In the right panel of Figure [3l we show the relative inten- 
sity profile of the NIR scattered light derived from the 
calculations with /enhance = 6.6 but /g a p(r) = 1 every- 
where (i.e., no gap). We compare the calculations with 
and without self-gravity, but there is no significant dif- 
ference between the two calculations. 

As long as the disk is massive enough and there exists 
a gap structure in the surface density, the bump struc- 
ture in the NIR profile is always present regardless of 
the location of the gap or the underlying surface density 
structure of the disk. Figure 0] shows the NIR profile for 
the model with r gap = 150AU. In this figure, we compare 
the results with /enhance = 1 and /enhance = 6.6. We see 
that there is a bump structure just inside (and also out- 
side) of the gap region. In Figure El we show the results 
of the NIR profile for the disk model with E a r" 1 . The 
disk model is the same with equations (UJ)-© except for 
the power of the surface density distribution, and we as- 
sume a = 0.95, r gap = 100AU and A gap = 30AU for the 
gap profile. We compare the results with /enhance = 1 
and /enhance = 6.6. In this model, again, we see a bump- 
like structure inside and outside of the gap region for 

/enhance — 6.6. 

We now discuss the origin of the bump structure in 
the NIR observation. In the NIR observation, we mainly 
observe the light from the central star scattered at the 
surface of the disk. We define the surface of the disk 
as the place where the optical depth towards the central 
star, T rac j, becomes unity. 

The strength of the scattered light, which is given by 
equation (|9]), can be well approximated by using the 
grazing angle, which is the angle between the incident 
ray from the central star and the dis k surface (see e.g., 
Ilnoue et al.ll200i IJang-CondelHl2009l) . The grazing an- 
gle can be calculated from the density profile. We have 
checked that the scattered light profile is indeed well re- 
produced by the grazing angle times a geometrical factor. 
Therefore, we discuss the effects of self-gravity on the be- 
havior of the surface of the disk and the grazing angle in 
the following paragraphs. We note that all the scattered 
light profiles given so far are calculated by directly inte- 
grating equation (j9|). 

We first consider the case where the self-gravity of the 
disk is negligible. In Appendix, we show explicit calcu- 
lations of the disk surface and the grazing angle in the 
case where vertical hydrostatic equilibrium is reached in 
a non-self- gravitating disk. It is shown that (1) the more 
the disk material, the high er in the altitude the disk sur- 
face is (see equation (|A7[| ). and (2) grazing angle ft is 
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approximately given by 

P~h-pT, (10) 

where h = c/ril is the aspect ratio of the disk, p Q is 
the density at the disk midplane, and (po) is the radial 
average of the midplane density calculated by 

<p ) = \ J dr' Po (r') (11) 

(see equation (|A12|) V If there exists a gap in the ra- 
dial surface density profile and the width of the gap is 
small compared to the disk scale, the midplane density 
at the gap location is low while the radially averaged 
density is not significantly affected by the presence of 
the gap. Therefore, we conclude that the scattered light 
traces the surface density profile well, thereby the NIR 
scattered light is well correlated with sub-mm dust con- 
tinuum emission. It is possible to see in the left panel of 
Figure [2] that the location of the dark region in the NIR 
profile and the location of the gap in the surface density 
profile are very close each other. 

If we include the effects of self-gravity, it acts to pull 
the disk material towards the disk midplane. Therefore, 
the disk self-gravity in general lowers the disk surface. 
Therefore, we expect that the position of the disk sur- 
face is maximized at a certain value of Q = Q c . Figure 
[6] shows the disk surface at 50AU (it is well away from 
the gap region) for various models with /enhance- The 
horizontal axis shows the disk Q-value at 50AU. It is 
shown that the disk surface is at the highest place when 
Q = Qc ~ 5. The exact value of Q c depends on the place 
of the disk we consider. For example, we have checked 
that the value of Q c is approximately 3 at 200AU. How- 
ever, the overall behavior of the location of the disk sur- 
face as a function of the disk mass is similar. 

Let us now consider what Figure |5] indicates for one 
particular disk model, in which there is a rapid variation 
of Q-parameter in the radial direction. From Figure [H 
it is possible to say that if there exists a place where 
the disk Q-parameter changes rapidly and crosses Q ~ 
Q c , we expect that the disk surface is at the highest 
place at such radius. The grazing angle is maximized 
there, and therefore, the scattered light from such radius 
is prominent. Figure [7] shows a schematic picture of this. 
In our disk model, there is a significant change in the 
disk Q-value at the location of the gap, and at the edge 
of the gap, there exists a place where Q ~ Q c . Figure |5] 
compares the profile of the grazing angle in the case of 
/enhance = 1 (almost non-self-gravitating) and /enhance = 
6.6 (self-gravitating disk). It is clearly shown that the 
grazing angle is maximized at the gap inner edge in the 
case where self-gravity becomes important. 

The appearance of the maximum of the grazing angle 
at the edge of the gap depends on how massive the disk 
is and to what extent self-gravity effects are important 
in the vertical structure of the disk. The bump is more 
enhanced if the disk is more massive. In order to see 
this, we calculate the strength of the peak at the inner 
edge of the gap as follows. First, we fit the brightness 
profile between 35AU and 50AU by power law. We have 
chosen this region from our data since this region is not 
significantly affected by the presence of the gap (see Fig- 



ure 13]). Then, we normalize the profile with the power 
law profile obtained. Finally, we look for the maximum 
between 60AU and 100AU from the normalized profile 
to find the location and the strength of the peak. The 
choice of 60AU and 100AU is arbitrary, but the point 
here is that the bump in NIR profile is included in this 
region. 

Figure [9] shows the relation between the peak strength 
and the averaged Q-value between 35AU and 50AU. The 
models with various values of /enhance and dust extinc- 
tion are shown. The different points sharing the same 
Q-value correspond to the different values of dust extinc- 
tion at NIR. It is shown that the smaller the Q-value, the 
stronger the peak, and that there is only small range of 
the Q-value that can cause the peak strength of 1.15 — 1.2. 
There is no significant peak for low-mass disks. For the 
particular model of /enhance = 6.6, which we have dis- 
cussed in detail in the previous section, the peak strength 
is of the order of 10%. 

In summary, the bump structure at the edge of the 
gap in the NIR scattered light profile is a characteris- 
tic feature arising from the effects of the self-gravity. It 
is also noted that, in the discussion above, the essential 
assumption is the vertical hydrostatic equilibrium, rapid 
variation of the Q-parameter, and the dominance of the 
scattered light in NIR observations. We therefore expect 
that the profile of the surface density is not necessarily a 
"gap" in order for the bump in the NIR imaging observa- 
tions to be produced. It is not necessarily axisymmetric, 
and it can be a "hole" or other such features that exhibit 
the change in the Q-parameters, as long as the hydro- 
static equilibrium is reached. 

4. DISCUSSION: NIR IMAGING OBSERVATIONS OF 
HIGH-MASS DISKS 

We have seen that the bump structure at the gap edge 
in the NIR observations reflects the fact that the disk 
self-gravity comes into play in determining the vertical 
structure of the disk. We now discuss the implications of 
the results to the observations. 

We first note that the peak strength shown in Figure [9] 
can, in principle, be calculated from the observed relative 
intensity profile only. From Figure O it may be possible 
to say that if we can detect the peak at the edge of the 
gap, it may be possible to constrain the disk Q- value 
more effectively than by using the sub-mm continuum 
emission profile alone, which always has an uncertainty 
of at least a factor of two arising from the uncertainty of 
the dust parameters. 

It is noted, however, that this argument would apply 
for the disk with a limited range of mass. For low-mass 
disks (having high Q-value), the peak is not very signif- 
icant and it is difficult to constrain the disk mass. This 
method of estimating the disk mass can be useful for the 
disk with the Q- value less than 4 or around. It is also 
noted that for very high-mass disks with Q-value less 
than 2 or around, the disk turbulence is important and 
the vertical hydrostatic balance may not be reached. 

In order to see this self-gravity effects in the actual ob- 
servations, it is necessary to have a very good sensitiv- 
ity in the NIR direct imaging observations. As we have 
seen in Figure |H1 it is necessary to detect the peak with 
the strength of 5 to 10 percent when the disk Q-value 
is about 3-4. The typical error of the NIR disk imag- 
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ing observations is rather difficult to estimate, since it 
depends on how well the stellar light can be suppressed. 
In the specific example of the obs ervation of the disk 
around AB Aur by AEOS telescop e (iQppenheimer et"al~l 
l200Sh . Uang-Condell and Kuchnen {2010) estimate that 
the error is typically at the level of ~ 20% by fitting the 
intrinsic azimuthal variation by a certain fitting function 
and by calculating the rms deviation. It is above the level 
of the peak strength calculated in this paper. However, 
in future observations, better suppression of the speckle 
noise may be made possible. 

For the particular model presented in this paper, the 
required resolution of the telescope can be reached by the 
current instrument like Subaru. We have checked this 
by assuming the disk is at 140pc and convolve the NIR 
profile with the Gaussian with FWHM 0.04 , which is the 
typical resolution of the Subaru telescope. The resulting 
image is not much different from the model image. 0.04 
at 140pc corresponds to 5.6AU, but the model image 
varies with the spatial scale of at least 10AU. 

We note that although we show that there may be a 
feature unique to a self-gravitating disk in the NIR scat- 
tered light observations, it is essential to compare this 
feature with the sub-mm observations. Since we can ob- 
serve only the tenuous surface of the disk in the NIR 
observations, small surface structure, such as a blob, can 
also cause the bump-like structure which ca sts a shadow 
over the disk ( Jang-Condcll and Bossl l2007l) . It is there- 
fore necessary to have a sub-mm observation which sug- 
gests that the disk is truly massive and there exists a gap 
in the surface density. The correlation between NIR and 
sub-mm images must also be checked in order to assure 
that the hydrostatic balance in the disk vertical direction 
is reached. 

Sub-mm observations are also important in determin- 
ing how the bump is produced. Another mechanism 
that may cause the bump-like feature in the NIR scat- 
tered light may be localized pressure bumps or a mi- 
grating protoplanet. A pressure bump may appear, 
for example, as a result of b aroclinic instability (e.g., 
iKlahr and Boswnheimerll2003l ). A protoplanet migrating 
inward in a disk may accumulate gas at the inner edge 
of a g ap and therefo re, a bump-like feature may appear 
(e.g., iRafikovl [20021) . These processes both act even if 
the disk mass is low enough and the effects of self-gravity 
is negligible. Since both processes involve the accumula- 
tion of the gaseous materials, such bump structures may 
also appear in the sub-mm dust continuum emission. It 
is therefore stressed again that the combination of the 
sub-mm and the NIR observations is important. 

The sub-mm imaging observations should also be used 
to make a disk model that can be compared with NIR 
observations. We note that the results presented in Fig- 
ure [9] are based on the specific model used in this paper. 
Although we expect that the general trend of the rela- 
tion between the Q-value and the peak strength is not 
very much different if we incorporate other forms of the 
gap profile or other disk models, it is necessary to make 
a model that fits with the sub-mm observation data of 
the surface density profile. It is essential to have at least 
the same level of spatial resolution between the NIR and 
sub-mm telescopes. 

In summary, if the disk mass is high and the Q- value is 



less than ~ 4, the effects of self-gravity would come into 
play in determining the disk vertical structure in hydro- 
static equilibrium, and this effect would appear in future 
NIR direct imaging observations as a detailed structure 
such as a bump and a gap, if we have sufficient sensitivity. 
The interpretation of the NIR imaging data would there- 
fore become complicated, especially in high-mass disk 
case. However, such detailed features may carry impor- 
tant information about the physical states of the disk 
such as disk mass. The combination between the sub- 
mm imaging observations and NIR imaging observations 
are also important in the interpretation of the image. 

5. SUMMARY AND CAVEATS 

In this paper, we have investigated the effects of self- 
gravity on the vertical structure of a protoplanetary disk 
with a gap, and discussed the possible signatures of a self- 
gravitating disk in direct imaging observations. We have 
seen that the self-gravity of the disk causes a bump fea- 
ture to appear at the edge of the gap in the NIR imaging 
observations, especially for the disks with Q-parameter 
with ~ 3 — 4. We have discussed that, although the 
interpretation of the NIR image would become compli- 
cated due to the effects of self-gravity if the disk is very 
massive, detailed analyses of the morphology of the ob- 
served image would provide us useful information about 
the physical state of the disk. 

The combination of the sub-mm and NIR images is 
important to assure that the disk is truly massive and is 
in hydrostatic equilibrium in the vertical direction, and 
to check whether other possible causes of the bump may 
play a role. Also, since the strength of the dust con- 
tinuum emission gives a rough estimate of the surface 
density, it is possible to check the consistency of the 
model (i.e., /enhance may be estimated both from the 
peak strength and the strength of the dust continuum 
emission) if such feature is detected. 

Although we have considered a specific disk model with 
an axisymmetric gap, we expect that the bump structure 
associated with the gap may be a universal feature in a 
self-gravitating disk that has a rapid variation in the Q- 
value. This feature comes from the fact that the surface 
of the disk is at the highest altitude in the vertical direc- 
tion at a certain value of Q-parameter. The ring geome- 
try is not essential in producing such a bump associated 
with a gap. 

We note that we have assumed the vertical hydrostatic 
equilibrium, and therefore, our model cannot be used for 
a turbulent disk. The turbulence induced by the self- 
gravity is expected to occur when the disk Q-value be- 
comes less than ~ 1.5 — 2, and in this case, the structure 
in the upper layer of the disk may be totally uncorrelated 
with that at the midplane. 

One of the main caveats of the model presented in this 
paper is the simplification of the temperature structure. 
Although the assumption of a vertically isothermal disk 
may not be seriously wrong for the outer disk, the irra- 
diation from the central star can cause a super-heated 
layer at the surface of the disk ( jChiang and Goldreichl 
Il997h . which can modify the observed signature in the 
NIR wavelength. We expect that the surface of the disk 
may be pumped up by the irradiation from the central 
star, and the bump structure we have discussed in this 
paper is more enhanced. Therefore, we expect that the 
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effects of irradiation may enhance the feature caused by 
the self-gravity. However, this issue is necessary to be 
investigated separately. 
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APPENDIX 
GRAZING ANGLE 

In this Appendix, we show the calculations of the grazing angle between the disk surface and the incident ray in 
the case of a non-self-gravitating disk. The objective of this section is to show explicitly that the observed NIR flux 
should be well correlated with the surface density of the disk, which is proportional to the observed sub-mm flux, if 
the following two conditions are satisfied: (1) the self-gravity of the disk is negligible and (2) the vertical hydrostatic 
balance is reached. 

If the hydrostatic equilibrium in the vertical direction is reached, the density profile of the non-self-gravitating, 
vertically isothermal disk is given by 

„2 



p(r,z) = p (»exp 



2H 2 



(Al) 



where H is the pressure scale height and po(r) is the density at the midplane. The midplane density is related to the 
surface density by 

E(r) 



Po(r) 



V2H{r) ' 



(A2) 



We first calculate the optical depth, r ra( j, between the central star and a point (r, z) in the disk. Since r ra( j is calculated 
by 

(A3) 

where / is the path connecting between the point (r, z) and the central star, we have 

,2 



Trad(r,z) = X 



\J r 2 + z 2 



dr'po(r') exp 



1 



2h(r') 2 r 2 



(A4) 



where h = H/r and we have assumed that the central star is a point source situated at the origin. 

We first look at the approximate position of the disk surface. We define the disk surface as a place where r ra( j = r s = 1. 
We assume that z <C r, which is justified later. The estimate of the order of magnitude of the integration in equation 
"TJ) yields 



dr'po(r') exp 



1 



and therefore, T rac i is 

Trad ~ XPore 

Therefore, the surface z s of the disk may be calculated by 



2h{r') 2 r 2 

-z 2 /2H 2 



p re 



2 /2H 2 



Hi / 2 log 



( XPor 



#4/2 log 



T S y/2h)' 



(A5) 
(A6) 

(A7) 



If we use So = 10g/cm 2 , X = 100cm 2 /g, t s = 1 and h = 0.1, we have z s ~ 4H. Therefore, the disk surface resides 
at several scale heights above the disk midplane. Since H <C r in a protoplanetary disk, this confirms the initial 
assumption of z <C r. 

We now look at the grazing angle ft between the disk surface and the incident ray. Since z s <^ r, f3 can be calculated, 
to a good approximation, by 



* 
ar 



surf 



r 



where dz/dr should be evaluated along the disk surface, 

dz 
dr 

surf 



dr Iad /dz 



(A8) 



(A9) 
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It is to be noted that from equation (|A4[) . T ra d(r, z ) lS & function of r and z/r only. The grazing angle j3 is calculated 

as 



yjl+(z/r) 2 p (r) exp[(l/2/i 2 )(z/r) 



d u [vTT^ 2 f dr' po(r-')exp[-u 2 /2/i 2 ]] 
Since z s /r ~ /i <C 1, the denominator can be approximated as 



s/r 



1 - ± ^P 1 / r ^Vo(r')%^e-« 3 / 2 MO 



Therefore, the grazing angle is given by 

z. 



h(r' 



p (r)exp[-(l/2/ l 2 )(z s /r) 2 



r r- 1 J r dr'p (r')(z s /rh(r')) 2 exp[-l/2/i(r') 2 (z s /r) 2 ] ' 
Since z s /r ~ h, the order of the magnitude of /3 may be 

Po 



(1/r) / dr'p (r')' 



(A10) 



(All) 



(A12) 



(A13) 



Therefore, we conclude that if there exists a gap in the disk and the gap width is narrow compared to the disk scale, 
the value of /3 is almost proportional to the midplane density, which is proportional to the surface density since we 
assume the hydrostatic equilibrium in a non-self-gravitating disk. Therefore, in this case, the NIR scattered light, 
whose intensity is proportional to the grazing angle, should be well-correlated with the sub-mm thermal emission, 
which reflects the profile of surface density. 
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Fig. 1. — An example of the surface density profile and the disk Q-parameter profile in the model. The left panel shows the surface 
density of the gas with /enhance = 1 in equation |TJ, and the right panel shows the disk Q-parameter for this model. 
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Fig. 2. — The radial profile of NIR observation (left) and sub-mm (right) observation. The surface brightness normalized by that at 
50AU is shown. The solid line corresponds to the disk model with /enhance = 1 (low-mass model) and the dashed line is for /enhance = 6.6 
(high-mass model). The results with extinction \ = 100cm 2 /g (NIR) and \ = 0.003cm 2 /g (sub-mm) is shown. In the left panel, we also 



show the profile of the surface density in the arbitrary unit by the dot-dashed line, 
mainly discussed in this paper. 



The arrow indicates the "bump" structure which is 
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Fig. 3. — Comparison between the calculations with self-gravity and without self-gravity. The left panel shows the disk model with a 
gap, and the right panel shows the disk model without a gap. The solid line shows the model with self-gravity, and the dashed line is the 
model without self-gravity. We use the model with /enhance = 6.6, which corresponds to the "high-mass model" in Figure [2] 
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Fig. 4. — The results of NIR profile for the model with the gap at 150AU. The models with /enhance = 1 (solid line) and with 
/enhance = 6.6 (dashed line) arc shown. The bump inside and outside the gap region is seen in the case of the high mass disk. 
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Fig. 5. — The results of NIR profile for the model with the surface density profile of E oc r . The models with /enhance = 1 (solid line) 
and with /enhance = 6.6 (dashed line) are shown. The bump inside and the outside the gap region is seen in the case of the high mass disk. 
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Fig. 7. — Schematic diagram for the origin of bump at the gap edge in NIR observation. In the region where the Q-value is small, the 
surface is a relatively lower place because of the effects of the self-gravity. In the region of the gap, the surface is also at a lower place 
because of the lack of the material. For the intermediate values of Q, the surface is at the highest place and therefore, the grazing angle is 
maximum, which leads to the formation of bump in the scattered light. 
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Fig. 8. — The grazing angle (in radian) as a function of radial distance (in AU) in two different / cn hancc parameters. The solid line 
shows the case of / cn hancc = 1-0 and the disk self-gravity is almost negligible. The dashed line shows the case of /enhance = 6.6. The effect 



of self-gravity appears at the gap edge, and there is a maximum of the grazing angle around 70AU region. The results with k : 
are shown. 
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Fig. 9. — The relation between the strength of the peak of the bump in the NIR observations and the Q-value. The Q-values are averaged 
over the radius between 35 AU and 50 AU. Different points with the same values of Q correspond to different extinction coefficients \ m 
the NIR wavelength. 



